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We extend the concept of graph isomorphisms to multilayer networks with any number of “as¬ 
pects” (i.e., types of layering). In developing this generalization, we identify multiple types of iso¬ 
morphisms. For example, in multilayer networks with a single aspect, permuting vertex labels, layer 
labels, and both vertex labels and layer labels each yield different isomorphism relations between 
multilayer networks. Multilayer network isomorphisms lead naturally to defining isomorphisms in 
any of the numerous types of networks that can be represented as a multilayer network, and we 
thereby obtain isomorphisms for multiplex networks, temporal networks, networks with both of 
these features, and more. We reduce each of the multilayer network isomorphism problems to a 
graph isomorphism problem, where the size of the graph isomorphism problem grows linearly with 
the size of the multilayer network isomorphism problem. One can thus use software that has been 
developed to solve graph isomorphism problems as a practical means for solving multilayer network 
isomorphism problems. Our theory lays a foundation for extending many network analysis methods 
— including motifs, graphlets, structural roles, and network alignment — to any multilayer network. 
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I. INTRODUCTION 

Network science has been very successful in investigations of a wide variety of applications in a diverse set of 
disciplines. In many situations, it is insightful to use a naive representation of a complex system as a simple, binary 
graph, which allows one to use the powerful methods and concepts from graph theory and linear algebra; and numerous 
advances have resulted from this perspective [1]. As network science has matured and as ever more complicated data 
have become available, it has become increasingly important to develop tools to analyze more complicated graphical 
structures [2, 3]. For example, many systems that were typically studied initially as ordinary, time-independent 
graphs are now often represented as time-dependent networks [ ], networks with multiple types of connections [5] , or 
interdependent networks [6] . Recently, a multilayer-network framework was developed to represent a large number of 
such networked systems [2], and the study of multilayer networks has rapidly become arguably the most prominent 
area of network science. It has achieved important results in a diverse set of fields, including disease dynamics [7], 
functional neuroscience [8], ecology [9], international relations [10], transportation [11], and more. 

With the additional freedom in representing a multilayer network, numerous ways to generalize network concepts 
have emerged [2, 3]. The different definitions can arise from different modeling choices and assumptions, which also 
have often been implicit (rather than explicit) in many publications. To make sense of the multitude of terminology 
and develop systematic methods for studying multilayer networks, one needs to start from first principles and define 
the fundamental concepts that underlie the various methods and techniques from network analysis that one seeks 
to generalize. For example, exploring the fundamental question, “How is a walk defined in multilayer networks?”, 
led to breakthroughs in generalizing concepts such as clustering coefficients [12, 13], centrality measures [14-16], and 
community structure [17-19] in multilayer networks. In this article, we answer another fundamental question: “When 
are two multilayer networks equivalent structurally?” by generalizing the concept of graph isomorphism to multilayer 
networks. 

Any attempt to generalize a method that relies on graph isomorphisms to multilayer networks also necessitates 
generalizing the concept of graph isomorphisms. Very recently, there has been work on methods relying on (some 
times implicit) generalizations of graph isomorphism—especially in the context of small subgraphs known as “motifs” 
—for many network types that can be represented as multilayer networks [21-26]. Further, other network analysis 
tools, such as structural roles [27, 28] and network comparison methods [29-33], are based on graph isomorphisms. 

Defining isomorphisms for multilayer networks yields isomorphism relations for each of the wide variety of network 
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types that can be expressed using a multilayer-network framework. For example, one obtains isomorphisms for 
multiplex networks (in which edges are colored), interconnected networks (in which vertices are colored), and temporal 
networks [2], Instead of defining isomorphisms and related methods and tools separately for each type of network, 
we develop a general theory and set of tools that can be used for any types of multilayer network [34]. With our 
contribution, we hope to avoid a confusing situation in the literature in which elementary concepts, terminology, tools, 
and theory are developed independently for the various special types of multilayer networks. 

The rest of this paper is organized as follows. In Section II, we introduce the basics concepts, lay out the ideas 
behind multilayer isomorphisms, and summarize the results of our article. In Section III, we give the permutation- 
group formulation of multilayer network isomorphisms and enumerate some basic properties of multilayer network 
isomorphisms and related automorphism groups. In Section IV, we show how to solve a multilayer network isomor¬ 
phism problem computationally by reducing it to an isomorphism problem in a vertex-colored graph. This reduction 
allows one to use graph isomorphism software packages to solve the multilayer network isomorphism problem, and we 
use it to show that multilayer isomorphism problems are in the same computational complexity class as the graph 
isomorphism problem. We provide tools for producing the reductions as a part of a multilayer analysis software [34]. 
In Section V, we give examples of how one can use multilayer network isomorphisms for multiplex networks, temporal 
networks, and interconnected networks. Finally, in Section VI, we conclude and discuss future research directions. 


II. BASIC CONCEPTS AND SUMMARY OF RESULTS 
A. Multilayer Networks 

In recent years, there has been a growing interest in generalizing the concept of graphs in various ways to study 
graphical objects that are better suited for representing specific real-world systems. This has allowed increasingly 
realistic investigations of complex networked systems, but it has also introduced mathematical constructions, jargon, 
and methodology that are specific to research in each type of system. The rapid development of such jargon has been 
overwhelming, and it has sometimes led to confusion and inconsistencies in the literature [2], 

To unify the rapidly exploding, disparate language (and disparate notation) and to bring together the multiple 
concepts of generalized networks that include layered graphical structures, the concept of a “multilayer network” was 
developed recently [2, 1 ]. Reference [2] includes a list of about 40 mathematical constructions that can be represented 
using the framework of multilayer networks. Most of these structures are variations either of graphs in which vertices 
are “colored” (i.e., vertex-colored graphs, see IV A and VB) or of graphs in which edges are “colored” ( i.e ., multiplex 
networks, see VA). Both types of coloring can also occur in the same system, various types of temporal networks 
admit a natural representation as a multilayer network [2], and other types of complications can also arise. This new 
unified framework has opened the door for the development of very general, versatile network concepts and methods, 
and the study of multilayer networks has rapidly become arguably the most prominent area of network science. See 
[2, 3] for reviews of progress in the study of multilayer networks. 



FIG. 1: Example of a multilayer network with two aspects. In this graphical structure, each entity can have an associated 
vertex-layer tuple in one or more layers, which are also organized into combinations of elementary layers. In this example, each 
layer has either A or B as the first elementary layer, and it has either X or Y as the second elementary layer. There are thus 4 
layers in total, and entities can have an associated vertex-layer tuple in one, two, three, or all four layers. [This plot is inspired 
by a figure in [2].] 

The formal definition of a multilayer network needs to be able to include the various layered network structures in the 
literature. In Ref. [2], we (and our collaborators) defined a multilayer network as a quadruplet M = ( Vm,Em , V, L) [2]: 
The set V consists of the vertices of a network, just as is an ordinary graph. Each vertex resides in one or more 
uniquely-named layers that are combinations of exactly d elementary layers , where each of these elementary layers 
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corresponds to an “aspect”. That is, each aspect is a different type of layering. For example, a social network that 
changes in time and includes social interactions over multiple communication channels has two aspects—one for time 
and the other for the type of social interaction -and so a layer represents one type of social interaction at a given 
time. The sequence L = {L a }f =1 consists of the sets of elementary layers for each of the d, aspects, and we use the 
symbol L = L\ x • • • x Ld to denote the set of layers. Each vertex can be either present or absent in a layer, and 
we indicate the presence of a vertex by including its combination with the layer in the set Vm C f x L\ x • • • x Ld 
of vertex-layer tuples. Finally, we define the set Em C Vm x Vm of edges between pairs of vertex-layer tuples as in 
ordinary graphs. See Fig. 1 for an example of a multilayer network with two aspects. 


B. Isomorphisms in Graphs and Multilayer Networks 

A graph isomorphism formalizes the notion of two graphs having equivalent structures. The structure is what is 
left in a graph when one disregards vertex labels. That is, two graphs are isomorphic if one can transform one graph 
to the other by renaming the vertices in one of the graphs. Note that the edges do not have their own labels but they 
are determined by the vertex labels of the two endpoints, and those labels are also updated in the transformation. 

To be able to give a mathematical definition of a graph isomorphism, we first define a vertex map as a bijective 
function 7 : V —> V' that relabels each vertex of the graph G = (V, E) with another distinct label. We use the 
following notation to relabel vertices of G using 7 : 

(1) F 7 = { 7(w ) I v G H} ; 

(2) E 1 = {( 7 (f), 7(u)) I (v,u) € E}; 

(3) G 7 = (H 7 ,# 7 ). 

With this notation, two graphs G and G' are isomorphic if there exists 7 such that G 7 = G'. 

One can define isomorphisms for multilayer networks in very similar manner. The idea is again that two networks 
are equivalent structurally if the vertices in one of them can be relabeled so that the first network is turned into 
exactly the second one. To do this, we need some additional (and slightly more cumbersome) notation: 

(!) v m = {(7(u),a) I (v,a) e V M }; 

( 2 ) e m = {((7(«),a),(7(w)i^)) I ((v,a),(u,/3)) € E M }; 

(3) M 7 = (V^, Ejj, V 7 , L). 

Note that (v. a) = (v, aq,..., a^), where a is a vector of layers. With the above definitions, we can now say that two 
multilayer networks M and M' are vertex-isomorphic when there exists 7 such that M 7 = M'. 

A vertex isomorphism is a natural extension of the standard graph isomorphism to multilayer networks, but it is 
not the only one. In a vertex isomorphism, one disregards only the vertex labels (but retains the layer labels) when 
comparing two multilayer networks. This choice is justifiable in some applications, but in others one might wish to 
also disregard the layer labeling. For example, one can map temporal networks into multilayer networks so that each 
time instance is a layer [ 2 ], and in this case two temporal networks are vertex-isomorphic if ( 1 ) the network has the 
same structure and order of structural changes and (2) the exact timings the structural changes are equal. However, 
if one is interested only in the relative order of the changes that take place in the network, one needs to be able to also 
disregard the layer labels. To do this, one can proceed in very similar way as for vertices, as it requires a function to 
relabel the layers. Specifically, we say that a bijective function S a : L a —> L' a is an elementary-layer map that renames 
the elementary layers of a network. 

One also may want to be able to relabel all of the elementary layers or only a subset of them in a multi-aspect 
multilayer network. We define a function S : L —)• LI that relabels all elementary layers, 8(a.) = (5i(ai),..., 8d(ctd )), 
and call it a layer map. A partial layer map 8 3l is a layer map for which 5 a = 1 if a 7 ^ ji for all l and where 1 is 
an identity map. That is, a partial layer map only relabels elementary layers that use some subset of all aspects. We 
say that these aspects ji are “allowed” to be mapped. We are now ready to define notation that formalizes the above 
ideas of how layer maps affect multilayer networks: 

(1) 1/ = and = {«5 a (a) \ a e L a } ■ 

( 2 ) v m = {(<b <5( a )) I (v, a) e V M } ; 

(3) E S M = {((v,8(a)),(u,8(f3))) \ ((v, a), (u, (3)) £ V M } ; 
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FIG. 2: Four examples of multilayer networks with one aspect: (a) M a , (b) Mb, (c) M c , and (d) Md■ The multilayer network 
M a is vertex-isomorphic to Mb, because there is a permutation 7 = (12 3) of vertex labels such that M2 = Mb- We can 
thus write M a =0 Mb- The network M a is layer-isomorphic to M c , and we write M a =1 M c because there is a permutation 
5 = (X Y) of layer labels such that M^ = M c . The network M a is also vertex-layer isomorphic to Md, and we write M a = 0,1 Md 
because there is a vertex-layer permutation ( = ("1,3) such that Mf = Md- Note that M a is not vertex-isomorphic to M c or 
Md, and it is not layer-isomorphic to Mb or Md- However, M a is vertex-layer isomorphic both to Mb and to M c . 

(4) M s = {V^,E s m ,V, I/). 

We can now say that two multilayer networks are layer-isomorphic when there exists a 6 such that M s = M'. Because 
of the intrinsic complications in defining general multilayer networks with any arbitrary number of aspects, the above 
notation is a bit cumbersome. In Section III A, we will make the notation less cumbersome, at the cost of also making 
it less explicit. 

We have defined isomorphisms related to relabeling either vertices or layers, but there is no reason why one cannot 
simultaneously do both of these. We thus define the vertex-layer map C = ('y, <5i,..., 6d) as a combination of a vertex 
map 7 and a layer map S. A vertex-layer map acts on a multilayer network such that a vertex-map and layer-map act 
sequentially on the network: M = (M 7 )" 5 . Clearly, the order in which the vertices and layers are relabeled does not 
matter, and vertex maps and layer maps commute with each other, so ( M 7 ) <5 = (M* 5 ) 7 . The vertex-layer maps can be 
used to define vertex-layer isomorphisms in the same way as one defines vertex isomorphisms and layer isomorphisms. 

We now collect all of our definitions of multilayer-network isomorphisms. 

Definition II. 1. Two multilayer networks M and M' are 

(1) vertex-isomorphic if there is a vertex map 7 such that M 7 = M '; 

(2) layer-isomorphic if there is a layer map d such that M s = M '; 

(3) vertex-layer-isomorphic if there is a vertex-layer map £ = ( 7 , S) such that M^ = M'. 

Layer isomorphisms and vertex-layer isomorphisms are called partial isomorphisms if the associated layer maps are 
partial layer maps. 

We use the notation M =0 M' to indicate that networks M and M' are vertex-isomorphic. We indicate partial layer 
isomorphisms by listing the aspects that are allowed to be mapped (i.e., aspects that do not correspond to identity 
maps in the partial layer map) as subscripts. If the layer isomorphism is not partial, we list all of the aspects of the 
network. We use almost the same notation for vertex-layer isomorphisms, where the only difference is that we include 
0 as an additional subscript. For example, for a single-aspect multilayer network, =1 denotes a layer isomorphism 
and =0,1 denotes a vertex-layer isomorphism. For partial layer isomorphisms and vertex-layer isomorphisms, we use a 
comma-separated list in the subscript to indicate the aspects that one is allowed to map. For example, = 2 is a partial 
layer isomorphism on aspect 2 , and —0,1,3 signifies a vertex-layer isomorphism in which one is allowed to map aspects 
1 and 3 but for which the layers in aspect 2 (and in any aspects larger than 3) are not allowed to change. We will 
explain the reason for this notation in Section III A. 

We give examples of a vertex isomorphism, a layer isomorphism, and a vertex-layer isomorphism in Fig. 2. 


C. Summary of Results and Practical Implications 

1. Applications of Isomorphisms 

The idea of a graph isomorphism is one of the central concepts in graph theory and network science, and it is an 
important underlying concept for many methods of network analysis—including motifs [20], graphlets [31, 35], graph 
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matching [29, 30], network comparisons [31-33], and structural roles [27]. Defining multilayer-network isomorphisms 
thus builds a foundation for future work by allowing generalization of all of these ideas for multilayer networks. 
Multilayer-network isomorphisms can also be used to define methods and concepts that are not intrinsic to graphs. For 
example, one can classify multilayer-network diagnostics and methods based on the types of multilayer isomorphisms 
under which they are invariant. 


2. Applications of Automorphisms 

The modeling flexibility added by layered structures in networks has led to the discovery of qualitatively new 
phenomena (e.g., novel types of phase transitions) for processes such as disease spread and percolation [2, 36, 37]. It 
is interesting to examine how multilayer network architectures affect structural features such as the graph symmetries. 
One can study symmetries using automorphism groups of graphs, as these enumerate the ways in which vertices can 
be relabeled without changing a graph. We formulate the idea of graph automorphism groups for multilayer networks 
in Section III, and we introduce a simplifying notation in which we think of vertices as a “0th aspect”. We show 
that combining maps of different aspects preserves all of the symmetries that are present in these aspects, but that 
completely new symmetries can result from combining these maps (see Proposition III. 1) . For example, if a symmetry 
exists under vertex-isomorphism or layer isomorphism, it must also exist under vertex-layer isomorphism. However, 
a vertex-layer isomorphism can lead to symmetries that are not present under either vertex isomorphism or layer 
isomorphism. 

The multilayer network automorphisms that we define in the present work generalize notions of structural equiva¬ 
lence of vertices (or, more precisely, “role equivalence”, “role coloring”, or “role assignment”) [27, 28]. Other related 
notions of structural equivalences have been defined in specific types of multilayer networks. For example, in so¬ 
cial networks with multiple types of relations between vertices, one can study the “block models” that one obtains 
by considering different types of homomorpliisms [38-41]. Additionally, in coupled-cell networks (which can have 
multiple types of edges and vertices), the automorphism groups and groupoids—which one obtains by relaxing the 
global condition for automorphisms—have a strong influence on the qualitative behavior of dynamical systems on 
such networks [42-44]. 


3. Aspect Permutations 

In our definition of multilayer networks, the elementary layers are ordered, and it is important to note that this 
is simply for bookkeeping purposes. Additionally, one can think of the vertices as elementary layers of a 0th aspect: 
from a structural point of view, the vertices are the same as other types of elementary layers. This is evident from 
the definition of multilayer networks, but it is far from evident in typical illustrations, in which vertices and layers are 
visualized, respectively, as points and planes. One can permute the order in which elementary layers are introduced, 
and isomorphism relations remain the same as long as the aspects in which the renamings are allowed are permuted 
accordingly (see Section III C) . 


4- Practical Computations 

For practical uses, it is important that the various types of multilayer isomorphisms can be computed in a simple 
and efficient way. It is a standard practice to solve this type of computational problem by reducing the problem to an 
isomorphism problem in (colored) graphs by constructing auxiliary graphs and then applying existing software tools 
for finding graph isomorphisms [45]. The auxiliary graphs can become complicated as the number d of aspects grows, 
and a slightly different auxiliary graph construction procedure needs to be defined for all 2 d types of isomorphisms. 
In Section IV, we show how to construct such auxiliary graphs for general multilayer networks in a way that the size 
of the problem grows only linearly with the size of the multilayer network. This opens up a very straight forward and 
efficient way to apply our approach for practical data analysis of any kind of multilayer networks without requiring 
knowledge of reductions or explicit construction of auxiliary graphs. 


5. Application to Specific Network Types 

Most studies of multilayer networks usually consider specific types of multilayer networks rather than studying them 
in their most general form [2]. In Section V, we show how the theory of multilayer isomorphisms can be applied to 
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some of the most typical types of networks: multiplex networks, vertex-colored networks (i.e., networks of networks), 
and temporal networks. We also illustrate how the different implicit isomorphism definitions for temporal networks 
from the literature [21] are related to our multilayer isomorphisms (see Section V C). Our isomorphism definitions for 
multilayer networks are explicit, and anyone who is familiar with multilayer isomorphism can very easily transfer that 
knowledge to isomoprhisms in temporal networks. 

One of the most prominent use of graph isomorphisms is motif analysis, in which all subgraphs of a network 
are grouped into isomorphism classes and the numbers of subgraphs in each class are examined [20]. For both 
computational tractability and the ability to interpret the results of such calculations, such analysis typically relies 
on using a reasonably small number of isomorphism classes. This limits the sizes of subgraphs that are studied, as 
the number of isomorphism classes grows very rapidly as a function of number of vertices. Similarly, the number 
multilayer-network isomorphism classes grows very rapidly both as a function of the number of vertices and as a 
function of the number of layers. Consequently, the same limitations of motif analysis that apply to ordinary graphs 
also apply for multilayer networks. In Section V A, we examine the growth of the number of isomorphism classes in 
multiplex networks. This illustrates the type of compromise that one needs to make in the number of vertices and 
layers that can be considered in a subnetwork to ensure that the number of isomorphism classes is reasonable. 


III. PERMUTATION FORMULATION AND PROPERTIES OF MULTILAYER ISOMORPHISMS 

We now show how to formulate the multilayer-network isomorphism problem in terms of permutation groups, and 
we give some elementary results for multilayer-network isomorphisms and related automorphism groups. 


A. Permutation Formulation of Multilayer Isomorphisms 

We limit our attention (without loss of generality) to multilayer networks M in which each of the networks has the 
same set V of vertices and same sets {T a }^ =1 of elementary layers. [64] 

We can now formulate the isomorphism theory using permutation groups. Vertex maps are permutations acting 
on the vertex set V , and elementary layer maps are permutations acting on elementary layer sets L a . If we construe 
the group operation as the combination of two permutations, then all possible vertex maps form the symmetric group 
Sy, and all possible elementary layer maps for a given aspect a form another symmetric group Sl a (i.e., 7 £ Sy, and 
6a £ Slo)- The vertex-layer maps are given by a direct product of the symmetric groups of vertices and of elementary 
layers. 

For notational convenience, we define the set of vertices to be the “0th aspect” (i.e., we define L 0 = V). We 
also introduce the following notation for vertex-layer tuples: v = (v,a). By convention, we define subscripts for 
vertex-layer tuples so that vq = v and v a = a a for a > 0, where v £ V and a. £ L. It is also convenient to use 1 c 
to denote a group that consists of the identity permutation over elements of the set C. Additionally, recalling that 
C = ( 7 , 5i,..., Sd ), it is convenient to use the notation = £(v) = (j(v), S(a)) and v£ a = Ca(w)- 

We let p C {0,1,..., d} (with |p| > 1) denote the set of aspects that can be permuted. Given p, we can then define 
permutation groups 

P p = D p x---xD p , (1) 

where D? = S'i Q if a £ p and D? = t La if a <£ p. We denote the complementary set of aspects by p = {0,1 ,d}\p. 

We obtain vertex permutations forp = {0}, layer permutations when 0 ^ p, and vertex-layer permutations when 0 £ 
p and |p| > 1. Layer permutations or vertex-layer permutations are partial permutations if there exists a £ {1,... d} 
such that a ^ p. 

We can now define multilayer-network isomorphisms for a set of multilayer networks A4. 

Definition III.l. Given a nonempty set p, the multilayer networks M,M' £ A4 are p-isomorphic if there exists 
( £ P p such that M *» = M'. We write M = p M'. 

We denote the set of all isomorphic maps from M to M' by Iso p (M, M') = {£ £ P p : M = M'}. Similarly, we use 
Autp(M) = Iso p (M, M) to denote the automorphism group of the multilayer network M. 


B. Basic Properties of Automorphism Groups 

In Eq. (1), we constructed the groups P p as direct products of symmetric groups and groups that contain only an 
identity element. The automorphism groups are subgroups of these groups: Aut p (M) < P p . A permutation remains 
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in the automorphism group even if we allow more aspects to be permuted (i.e., if the set p is larger), and permutations 
that use only a given set of aspects are independent of permutations that use only other aspects. We formalize these 
insights in the following proposition. 

Proposition III.l. The following statements are true for all M and \p\ > 0: 

(1) Aut P '(M) < Aut P (M) if p' C p • 

(2) Aut Pl (M)Aut P2 (M) < Autp(M) if P\,P 2 Q P, with pi flp 2 = 0; 

(3) = 1 ==> C ^ = 1 f or a ll i if C W G Aut Pi (M) and pi C\pj = 0 for all i j . 

For a proof, see Section VIII A. 

It is important to observe in claim (2) of Proposition III.l that the subgroup relation can be proper even if p = piUp 2 . 
That is, the relationship Aut Pl (M)Aut P 2 (M) = Aut Pl u P 2 (Af) is not always true, but one can combine permutations in 
Pn and P P2 that are not in the automorphism groups Aut Pl or Aut P2 to obtain a permutation that is in Aut Pl Up2 (M). 
We give an example in Fig. 3. 



FIG. 3: An example demonstrating that one cannot always construct multilayer-network automorphism groups by combining 
smaller automorphism groups. That is, Aut pi (A/) Aut P2 (M) ^ Aut pl u P2 (A/) i n this example. In a directed multilayer network 
M with edge set Em = {[(1, X), (2, A)], [(1, X), (1, V)], [(2, V), (1, Y)], [(2, Y), (2, A)]}, both the vertex automorphism group 
Aut{ 0 }(lkf) and the layer automorphism group Aut{i}(M) are groups whose only permutation is the identity permutation, but 
the vertex-layer automorphism group A-itt{o,i}(M) has a permutation ((12), (XY)) in addition to the identity permutation. 


C. Aspect Permutations 

In the definition of multilayer networks, the order in which one introduces different types of elementary layers (i.e., 
aspects) only matters for bookkeeping purposes. For example, for a system that is represented as a multilayer network 
with two aspects, A and B , it does not matter if we assign index 1 to aspect A and index 2 to aspect B or index 1 to 
aspect B and index 2 to aspect A. The isomorphisms of type =i and = 2 in the former case become the isomorphisms 
of type =2 and =i in the latter case, and vice versa. Similar reasoning holds even if we consider the vertices to be a 
“0th aspect”, as we did in Section III A. 

To formalize the above idea, we introduce the idea of aspect permutations as permutations of indices of the as¬ 
pects (including the 0th aspect). We then show that multilayer-network isomorphisms are invariant under aspect 
permutations as long as the indices in the set p of aspects that are not restricted to identity maps are permuted 
accordingly. 

Definition III. 2. Let a € *S{o,...,d} be a permutation of aspect indices. We define an aspect permutation of a 
multilayer network as A a (M) = (Vf^, E' M , V 1 , L'), where 


o' 

1 

b 

II 

5? 

i-H 

Aa-i(d)) I v G V M } ; 

o' 

1 

b 

II 

( w <7- 1 (0)>---,' u <T- 1 (d))] 1 (v u ) e Em} 

(3) W = L ct - 1(0) ; 


(4) L' = 



See Fig. 4 for an example of an aspect permutation in a single-aspect multilayer network. For single-aspect multilayer 
networks, there is only one nontrivial aspect permutation operator, and we call the resulting multilayer network its 
aspect transpose. Multilayer networks that are vertex-aligned [2] (i.e., networks for which Vm = Vo X • • • x V,j) are 
often represented using adjacency tensors [2, 14, 46]. In this case, aspect permutations of multilayer networks become 
permutations of tensors indices [47-49] in the tensor representation. Note that aspect permutation is a meaningful 



(b) 


1 



FIG. 4: Two single-aspect multilayer networks, (a) M a and (b) Mb, that are aspect-transposes of each other: M a = A( 0 i)(Mf,) 
and Mb = 4l(oi )(M a ). The transposition operation preserves multilayer isomorphisms in the sense that a third multilayer 
network M c is vertex-isomorphic (respectively, layer-isomorphic) to M a if and only if A( 0 i) (M c ) is layer-isomorphic (respectively, 
vertex-isomorphic) to Mb. Similarly, M c is vertex-isomorphic (respectively, layer-isomorphic) to Mb if and only if A( 01 )(M C ) 
is layer-isomorphic (respectively, vertex-isomorphic) to M a \ and M c is vertex-layer-isomorphic to Mb (respectively, M a ) if and 
only if j4(oi)(M c ) is vertex-layer-isomorphic to M a (respectively, Mb). Illustrations produced using [34]. 


operation even for undirected multilayer networks, and it is different from the transpose operator, which reverses the 
orientations of the edges. 

Aspect permutations preserve the sets of isomorphisms as long as the indices in the isomorphism permutations are 
also permuted accordingly. 

Proposition III.2. The relation 

Iso P (M, M') = I q — i [Iso p-iMM), MM'))] (2) 

holds, where I a -i(C) = (Cr(o)> • • ■, C 7 (d)) an operation that permutes the order of elements in a tuple according to 

the permutation <r _1 and p a is a set in which each element of p is permuted according to the permutation a. 

For a proof, see Section VIII A. 


IV. SOLVING MULTILAYER ISOMORPHISM PROBLEMS 

To take full advantage of the theory of isomorphisms in multilayer networks, one needs efficient computational 
methods for finding isomorphisms between a pair of multilayer networks. One can proceed on a case-by-case basis 
for various types of networks, such as temporal networks [21], using standard techniques from the graph-isomorphism 
literature [45] . We will now use the same techniques to show how to reduce all of the multilayer-network isomorphism 
problems to vertex-colored-graph isomorphism problems. This reduction allows one to solve any kind of isomorphism 
problem for any type of multilayer network without the need to come up with and prove the correctness of a new 
reduction technique. 

In the reductions that we define, the size of the vertex-colored-graph isomorphism problem is a linear function of 
the size of the multilayer-network isomorphism problem and thus yields practical ways of solving multilayer-network 
isomorphism problems. We also use these reductions to show that solving multilayer-network isomorphism problems 
is in the same complexity class as ordinary graph isomorphism problems. This is unsurprising, as many generalized 
graph isomorphism problems are known to be equivalent [50], including ones that involve the very general relational 
structures defined in Ref. [51]. Another valid approach for our argument would be to reduce a multilayer-network 
isomorphism problem to other structures (e.g., to a fc-uniform hypergraph [52]), but the reduction to a vertex-colored 
graph yields practical benefits in terms of the ability to directly use software that is designed to solve isomorphism 
problems. 


A. Isomorphisms in Vertex-Colored Graphs 

A vertex-colored graph G c = (V c , E c , 7r, C) is an extension of a graph ( V c , E c ) with a surjective map n : V —> C that 
assigns a color to each vertex. We define a vertex map as a bijective map 7 : V c — S► V'. and introduce the following 
notation: Vff = {7(1;) | v £ V c }, Ef = {(7(1)), j(u)) | (v,u) £ E c }, 7t 7 (v) = 7r(7 _1 (u)), and Gf = (V 7 , Ef, 7r 7 , C). 
Two vertex-colored graphs, G c and G' c , are isomorphic if there is a vertex map 7 such that G 7 = G' c , and we then 
write G c = G' c . 
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For the purposes of isomorphisms, we can— without loss of generality—limit our attention to graphs with the 
vertex set V = {l,...,n}, where n is the number of vertices in the graph. This allows us to phrase the graph 
isomorphism problem in terms of permutations (similar to Section III A). The bijective map 7 in the definition of a 
graph isomorphism is again a permutation that acts on the set V of vertices, and the permutations form the symmetric 
group Sy- 

The vertex-colored-graph isomorphism problem is a well-studied computational problem, and several algorithms 
and accompanying software packages are available for solving it [45, 50, 51]. 

B. The Reduction 



FIG. 5: Example of a function f p that maps multilayer networks to vertex-colored graphs, (a) A multilayer network M 1 
with a single aspect, two layers, and three vertices, (b) The vertex-colored graph /{o}(Mi). One can use the mapping /{0} 
to find vertex isomorphisms in the multilayer network Mi. In other words, permutations of vertex labels are allowed, but 
permutations of layer labels are not allowed, (c) The vertex-colored network /{o,i}(Afi). One can use the mapping /{o,i} to 
find vertex-layer isomorphisms in the multilayer network Mi. In other words, both vertex labels and layer labels are both 
allowed to be permuted, (d) A multilayer network M 2 with two aspects, two layers in each aspect, and four vertices, (e) The 
vertex-colored graph f{i, 2 }(M 2 ). One can use the mapping /{i, 2 } to find layer isomorphisms in the multilayer network M2. In 
other words, permutations of layer labels are allowed in each aspect, but permutations of vertex labels are not allowed, (f) The 
vertex-colored graph /{!}(M2). One can use the mapping /{!} to find partial layer isomorphisms in the multilayer network 
M 2 . In other words, permutations of layer labels are allowed only in the first aspect, and permutations of vertex labels or layer 
labels are not allowed in the second aspect. 


The idea behind our reduction of multilayer-network isomorphism problems to the isomorphism problem in vertex- 
colored graphs is that we define an injective function f p such that two multilayer networks M and M' are isomorphic 
with a permutation from P p if and only if f p (M ) and f p (M') are isomorphic vertex-colored graphs. In this reduction, 
it is useful to consider the concept of an underlying graph Gm = (Vm,Em) of a multilayer network [2]. For two 
multilayer networks to be isomorphic, their underlying graphs need to be isomorphic. However, this is not a sufficient 
condition, because it allows (1) permutations in aspects that are not included in p and (2) permutations that occur 
in each layer independently of permutations that occur in other layers. Consider, for example, the multilayer network 
M a in Fig. 2 and the network M' that one obtains by swapping vertex labels 2 and 3 in layer X but not in layer Y. 
The underlying graphs Gm 0 and Gm 1 are then isomorphic even though there is no vertex-layer isomorphism between 
the two associated multilayer networks. 

We address the first issue above by coloring the vertices in the underlying graph so that its vertices, which correspond 
to vertex-layer tuples in the associated multilayer network, that are not allowed to be swapped are assigned different 
colors from ones that can be swapped. For example, for a vertex isomorphism in a single-aspect multilayer network, 
we color the vertices of the underlying graph according to the identity of their layers (i.e., by using a different color 
for each layer). We address the second issue above by gluing together vertex-layer tuples that share a vertex or an 
elementary layer by using auxiliary vertices. For example, for a vertex isomorphism in a single-aspect multilayer 
network, we add an auxiliary vertex in the underlying graph for each vertex v € V in the multilayer network, and 
we connect the auxiliary vertex to vertices in the underlying graph that correspond to v. This restricts the possible 
permutations: for each layer, one needs to permute the vertex labels in the same way. See Fig. 5 for an example of 
our reduction procedure. 

We define the reduction function f p for general M and p as follows. 

Definition IV. 1. We construct the reduction from multilayer networks to vertex-colored graphs f p : Ai —> Q c such 
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that f p {{V M , E m , V, L)) = (V G , E G , C, n) using 

(1) Vq = Vm U Vo, where the auxiliary vertex set Vo = |J agp > 

(2) E g = E m U E 0 , where E 0 = {(u Q , v)|v€ V M , a £ p} ; 

(3) C = pU Lp i x • • • x Lp m ; 

(4) 7 r(v g ) = a if v g £ L a and ir(v g ) = (v Pl ... % m ) if v g = v e V M ■ 

In addition to the reduction function f p that we need to solve the decision problem of two multilayer networks 
being isomorphic, we would like to be able to explicitly construct the permutations that we need to map a multilayer 
network to an isomorphic multilayer network. That is, we need a mapping between the permutations in multilayer 
networks and permutations in vertex-colored graphs. We define this map as follows. 

Definition IV. 2. Given a multilayer network M, we define the function g p from the permutations P p to permutations 
of vertex-colored graphs so that v g p ^ = v£ if v g £ Vm and v g p( ^ = v^ a if v g £ L a for any £ £ P p . 

The following theorem allows us to use f p and g p for the purpose of solving multilayer network isomorphism problems 
using an oracle for vertex-colored graph isomorphism. 

Theorem IV.l. Iso„(M, M') = g~ l [Iso(/ p (M), f p (M'))\ 

For a proof see Section VIIIB. 

From Theorem IV.l, it follows that one can also solve multilayer network isomorphism problems using the reduction 
to vertex-colored graphs that we have introduced. For example, one can use this reduction to determine if two 
multilayer networks are isomorphic, to define complete invariants for isomorphisms, and to calculate automorphism 
groups. We summarize these uses of Theorem IV.l in the following corollary. 

Corollary IV. 2. The following statements are true for all multilayer networks M, M' £ A4 and nonempty p: 

(1) M = p M’ ^ f P (M) “ f p (M'); 

(2) CoifpiM)) is complete invariant for = p if Cq is complete invariant for =; 

(3) Autp(M)=g-\Aut(f p (M))). 

For a proof, see Section VIIIC. 

We now define the “multilayer network isomorphism decision problem” and show that it is in the same complexity 
class with the graph isomorphism problem if one problem is allowed to be reduced to the other in polynomial time. 

Definition IV. 3. The multilayer network isomorphism problem ( MGI p ) gives a solution to the following decision 
problem: Given two multilayer networks M,M' £ At, is M = p M' true? 

The complexity class in which problems can be reduced to the graph isomorphism problem is denote here GI 
and many graph-related problems such as vertex-colored graph isomorphism problem and hypergraph isomorphism 
problem are known to be G/-complete [50]. 

Corollary IV. 3. MGI p is GI-complete for all nonempty p. 

For a proof, see Section VIIIC. 

We do the reduction from multilayer networks to vertex-colored graphs using the f p function that we defined 
earlier. We only need to show that this reduction is indeed linear (and thus also polynomial) in time. The reduction 
of graph isomorphism problems to multilayer network isomorphism problems is trivial if we allow the vertex labels to 
be permuted, because we can simply map the graph to a multilayer network with a single layer. If we cannot permute 
the vertex labels- i.e., if 0 ^ p —then we need to construct a multilayer network in which each vertex of the graph 
becomes a layer with only a single vertex and we then connect these layers according to the graph adjacencies. 
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V. ISOMORPHISMS INDUCED FOR OTHER TYPES OF NETWORKS 

In this section, we illustrate the use of multilayer network isomorphisms in network representations that can be 
mapped into the multilayer-network framework. As example, we use the three most common types of multilayer 
networks [2]: multiplex networks, vertex-colored networks, and temporal networks. In Section V A, we discuss iso¬ 
morphisms in multiplex networks. We focus on counting the number of nonisonrorphic multiplex networks of a given 
size (i.e., with a given number of vertices). In Section VB, we discuss isomorphisms in vertex-colored networks. In 
Section V C, we illustrate how multilayer network isomorphisms give a natural definition of the isomorphisms that 
are defined implicitly for temporal networks when analyzing motifs in them [21]. 


A. Multiplex Networks 
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FIG. 6: Isomorphism classes for multiplex networks with 3 vertices and 2 layers. We only include connected networks. We show 
vertex isomorphisms in the left panel and vertex-layer isomorphisms in the right panel. The count is the number of networks 
(with a fixed set of vertices and layers) that are mapped to each class. 

Multiplex networks have thus far been the most popular type of multilayer networks for analyzing empirical network 
data [2, 3]. One can represent systems that have several different types of interactions between its vertices as multiplex 
networks that are defined as a sequence of graphs {G a } a = {V a ,E a } a . It is almost always assumed that the set 
of vertices is the same in all of the layers V a = Vp for all a,/3 (although this is not a requirement), and multiplex 
networks that satisfy this condition are said to be “vertex-aligned” [2], 

One can map multiplex networks to multilayer networks with a single aspect by considering each of the graphs G a 
as an intra-layer network (i.e., a network in which the edges are placed inside of a single layer [2]). Optionally, one 
can add inter-layer edges (i.e., edges in which the two vertices are in different layers) by linking each vertex to its 
replicates in other layers. This is known as categorical coupling. Either using categorical coupling or leaving out all 
of the inter-layer edges leads to same isomorphism relations for multiplex networks. However, for ordinal coupling, 
in which only vertices in consecutive layers are adjacent to each other, the isomorphism classes can be different (see 
Section VC). A vertex isomorphism in multiplex networks allows the vertex labels to be permuted, but the types of 
edges are preserved. The layer isomorphism allows the types of edges to be permuted but only in a way that all of 
the edges of a particular type are mapped to a single other type. 

Analyzing small substructures using clustering coefficients in social networks and other multiplex networks have 
recently gained attention [12, 53-57]. Such structures have important (and fascinating) new features that go beyond 
clustering coefficients in ordinary graphs. Instead of there being only one type of triangle, there is very large number 
of different types of multiplex triangles and connected triplets of vertices. Such triadic structures have not been fully 
explored, though we discuss them in some detail in a recent paper [12]. Moreover, one can study larger subgraphs and 
induced subgraphs of multiplex networks by extending the analysis of “motifs” in graphs [20] to multiplex networks. 
There has already been interest in motif analysis in gene-interaction networks with multiple types of interactions 
[23], in food webs that can be represented using directed ordered networks [24], and in brain networks with both 
anatomical and functional connections [26]. Methods based on counting the number of isomorphic subgraphs, such 
as motif analysis, work best if the number of isomorphism classes is relatively small. Such analysis also necessitates 
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TABLE I: (Top) Number of isomorphism classes in multiplex networks for (left) vertex-layer isomorphisms and (right) vertex 
isomorphisms. (Bottom) The numbers of isomorphism classes with a given number of edges. All of the rows are symmetric 
around the maximum value(s), which we indicate in bold. The isomorphism classes were enumerated using [34]. 


the investigation of isomorphisms for their own sake, and they thereby serve as an important motivation (as well 
as an obvious future direction) for the present work. In Fig. 6, we enumerate all of the possible isomorphisms in 
connected multiplex networks with 3 vertices and 2 layers. We indicate each of the 16 vertex-isomorphism classes and 
10 vertex-layer-isomorphism classes. 

The problem of counting the nonisomorphic graphs that have some restrictions is known as the “graph enumeration 
problem” in graph theory, and such problems can be extended to multiplex networks (or multilayer networks in 
general) using the theory that we have introduced in the present paper. The number of undirected graphs with a 

fixed set of n vertices is 2W, and the number of nonisomorphic graphs also grows very quickly with n. In multiplex 
networks, the analogous problem is to count the number of multiplex networks with n vertices and b layers. For vertex- 
aligned multiplex networks, the number of networks is 2 b ( 2 ). In Table V A, we show the number of nonisomorphic 
vert ex-aligned multiplex networks for small values of n and b when considering vertex isomorphism or vertex-layer 
isomorphism. We produce the numbers in the table by systematically going through all of the networks of a certain 
size and categorizing them according to their isomorphism class [65]. The layer isomorphism problem for multiplex 
networks does not require one to solve the graph isomorphism problem, and it is easy to solve analytically. For layer 
isomorphisms, the number of nonisomorphic networks in a single-aspect vertex-aligned multiplex network is given by 

the formula ( 2 ^ 2 ^+ 6-1 )_ 


B. Vertex-Colored Networks 

One can represent networks with multiple types (i.e., colors, labels, etc.) of vertices using the vertex-colored graphs 
that we discussed in Section IV A. One can also map structures such as networks of networks, interconnected networks, 
and interdependent networks into the same class of multilayer networks [2], because one can mark each subnetworks 
in any of these structures using a given vertex color. 

One can map vertex-colored networks into multilayer networks by considering each color as a layer. One then adds 
vertices to the layer that corresponds to their color. Each vertex thus occurs in only a single layer, and one can add 
edges between the vertices in the resulting multilayer network exactly as they appear in the vertex-colored network. 
That is, in this multilayer-network representation, all inter-layer and intra-layer edges are possible. 

Vertex isomorphisms in this case are the normal isomorphisms of vertex-colored graphs, as vertex labels can be 
permuted but the colors are left unchanged. In layer isomorphisms, the vertex labels must be left untouched, but 
the colors can be permuted. For example, consider two networks with the same topology but different colorings 
that correspond to vertex classifications (e.g., community assignments [58]) of vertices. Two networks are then layer 
isomorphic if the two vertex classifications are the same. In a vertex-layer isomorphism, one can permute both the 
vertex names and the colors. 
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C. Temporal Networks 




FIG. 7: (a,b) Two event-based directed temporal networks that were used as an example in Ref. [21 ] to illustrate the difference 
between temporal motifs and flow motifs. The two temporal networks correspond to two distinct temporal motifs (i.e., two 
distinct isomorphism classes) but to a same flow motif (i.e., the same isomorphism class). The numbers next to the edges are 
times at which events take place. (c,d) Representations of the two temporal networks as vertex-aligned multiplex networks 
in which each vertex is present on each layer and the layers are ordinally coupled. This representation leads to the same 
isomorphism as used for temporal motifs in Ref. [21], and the two multilayer networks are not isomorphic, because the coupling 
edges fully determine the relative order of all layers. (e,f) Representation of the two temporal networks as non-vertex-aligned 
multiplex networks. In this representation, vertices are only present on layers in which they are active, and they are only 
adjacent to their replicas in other layers that participate in events. Consequently, similar to the isomorphisms that were used 
to define flow motifs in Ref. [21], the relative order of events is only important for events that are adjacent. The two multilayer 
networks constructed in this way are thus vertex-layer isomorphic. Multilayer-network illustrations produced using [34]. 


Temporal networks in which each edge and vertex are present only at certain time instances arise in a large variety of 
scientific disciplines (e.g., sociology, cell biology, ecology, communication, infrastructure, and more) [4], (One can also 
think about temporal networks with intervals of activity or with continuous time.) One can represent such temporal 
networks as multilayer networks [2, 14], although this is not the usual framework that has been used to study them. 
(See [17] for an early study that used this perspective.) Representing temporal networks as multilayer networks allows 
one to use ideas and methodology from the theory of multilayer networks to study them, and this has already been 
profitable in application areas such as political science [17], neuroscience [59], finance [60], and sociology [61]. More 
typically, one represents temporal networks either as contact sequences or time sequence of graphs [4]. Sequences of 
graphs are very similar construction to multiplex networks, where the key difference is that the order of the graphs in 
the sequence is important. One can map this type of temporal network to a multilayer network in very similar way as 
with multiplex networks. For temporal networks, however, one typically uses ordinal coupling instead of categorical 
coupling, although it is possible to be more general [2]. (In other words, instead of coupling all of the layer together, 
one only couples consecutive layers [2, 17].) 

A contact sequence consists of a set of triplets (u,v,t) that each represents a (possibly directed) contact between 
vertices u and v at time t. It is common to represent contact-sequence data as a sequence of very sparse graphs 
in which each distinct time stamp corresponds to a graph, and two vertices are adjacent in such a graph if they 
participate in an event at that time stamp [ ]. This representation leads naturally to the multiplex-like multilayer 
network representation of contact sequences that we described above. Alternatively, one can represent each event as 
a layer that only includes the two (or potentially more) vertices that participate in the event. The two vertices in the 
layer are each adjacent to its replicas in temporally adjacent layers. (See our earlier discussion of ordinal coupling.) 
These two alternative representations of temporal networks induce different isomorphism relations, and this difference 
is related to the difference between the temporal motifs and flow motifs from Ref. [21]. We illustrate this distinction 
using an example in Fig. 7. 

Contact sequences can also include delay or duration of the contact [4]. The delay (or latency) implies that the 
effect of a contact is not instantaneous. For example, in a temporal network of airline traffic, one can construe the 
flight time of each flight as a delay, and this can have an effect on the temporal paths and dynamical processes on the 
network [62] . This type of temporal network can also be represented using a multilayer-network framework [2] . For 
example, a flight that leaves city A at time t\ and arrives in city B at time £2 is represented as an edge from vertex 
A in layer t\ to vertex B in layer t2. Consequently, multilayer network isomorphisms can also be used for temporal 
networks with delays. 

In a network that is purely temporal, and which thus has only a single aspect, there are three different possible 
multilayer isomorphisms. (1) Two temporal networks are vertex-isomorphic if they exhibit the same temporal patterns 
at exactly the same time but between (possibly) different vertices. (2) Two temporal networks are layer-isomorphic 
if they exhibit exactly the same temporal patterns with exactly the same vertices, although the actual times (though 
not the relative order of events) can change. (3) Two temporal networks are vertex-layer isomorphic if they exhibit 












14 


exactly the same temporal pattern, though possibly using different vertices, but the actual times (although not the 
relative order of events) can change. 


VI. CONCLUSIONS AND DISCUSSION 

The theory of multilayer network isomorphisms illustrates the power of the multilayer-network formalism: Any 
concept or method that can be defined for general multilayer networks immediately yields the same concept or 
method for any type of network that can be construed as a type of multilayer network. The interpretation of the 
concepts or methods depends on the application and scientific question of interest, but the underlying mathematics is 
the same. In this sense, multilayer networks allow one to return to the early days of network science in which simple 
graphs were used to represent myriad types of systems and the same tools could be applied to all of them. The key 
difference is that multilayer networks allow one to represent much richer and application-specific structures. 

Going from graphs to multilayer networks adds a “degree of freedom” to ordinary networks (or multiple degrees of 
freedom if the number of aspects is larger than 1), and generalizing concepts defined for graphs thus typically leads 
to multiple alternative definitions [2]. This is also true for graph isomorphisms and any isomorphism-based methods 
in multilayer networks, and this is underscores why it is important to identify multiple types of multilayer network 
isomorphisms. Given a problem under study, one still needs to decide which of these generalizations to use. Naturally, 
one can also examine multiple types of isomorphisms. 

Our work on multilayer network isomorphisms lays the foundation for many future research directions in the study 
of multilayer networks. Motif analysis can now be generalized for any type of multilayer network once one defines 
a proper null model for the type of multilayer network under study. A good selection of network models already 
exist both for multiplex networks and for vertex-colored networks and similar structures [2] . Another straightforward 
application of isomorphisms in multilayer networks is the calculation of structural roles [27, 63] by defining two vertices 
to be structurally equivalent if they are equivalent under an automorphism. One can also examine other types of role 
equivalence in a multilayer setting. 

One of the challenges in isomorphism-based analysis methods is that they are computationally challenging even for 
ordinary graphs. We introduced a computationally efficient way of deducing if two multilayer networks are isomorphic 
and calculating multilayer network certificates by reducing the problem to the isomorphism problem for vertex-colored 
graphs. Although this method is efficient for general multilayer networks, there is room for improvement when one is 
only considering a specific type of multilayer network (such as multiplex networks). 

Our theory also forms a basis for methods that still need some additional work to be generalized for multilayer 
networks. For example, in interesting direction would be to define “approximate isomorphisms” or inexact graph 
matching [29] along with a way to measure how close one is to achieving an isomorphism. This would, in turn, allow 
one to define similarity measures between multilayer networks and techniques for “aligning” two multilayer networks. 
It would also make it possible to relax the conditions in role equivalence to better study structural roles in multilayer 
networks. 

Perhaps the most exciting direction in research on multilayer networks is the development of methods and models 
that are not direct generalizations of any of the traditional methods and models for ordinary graphs [2] . The fact that 
there are multiple types of isomorphisms opens up the possibility to help develop such methodology by comparing 
different types of isomorphism classes. We also believe that there will be an increasing need for the study of networks 
that have multiple aspects (e.g., both time-dependence and multiplexity), and our isomorphism framework is ready 
to be used for such networks. 
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VIII. PROOFS 

A. Proofs of basic properties of isomorphism and automorphism groups 

Proof of Proposition III.l. (1) Take any <f £ Aut p '(M). It follows that = M and £ £ P p because C £ P p >■ That 
is, C £ Autp(M). 

(2) Both Aut Pl (M) and Aut P2 (M) are subgroups of Aut p (M) because of (1). Their direct product is a group if 
they commute. Take any Q £ Aut Pl (M) and C,' £ Aut P2 (M). We have (CC')a = (Cl)a = (lC)a = (C'Qa if a £ p\, 
(CC')a = (lC')a = (C'l)a = (C'Oa if a £ P2, and (CC')o = (l)a = (C'Oa if a £ P1,P2- Therefore, CC' = C'C and 
Aut Pl (Af)Aut P2 (M) = Aut P2 (M)Aut Pl (M). 

(3) Let us look at arbitrary aspect a. Because pi n Pj = 0 for all i j, it follows that a is either a member of 

exactly one pi or of none of them. If a is not in any pi , then Dff = t^ a and fa ' 1 = 1 L a for all i. However, if a G Pj 
(i.e., the aspect a is in exactly one set), then (IliC^)a = (C (:d )a, and it thus follows that (C ( ^)a = 1 L a - Because a 
is arbitrary, we have shown that <f ^ = 1 for all i. □ 

Proof of Proposition III. 2. We first show that Iso P (M,M') C I a ~i [Iso p ^ (A a (M), A a (M’))]. We consider any £ G 
Iso p (M,M') and show that I a (C) £ Iso p a(A a (M), A a (M')). By a direct calculation, A a (M) Ia ^ = A a (M'): for 

vertex-layer tuples, (/ ct (Vm)) 7 " (c) = (o) > ’ ■ • ■> v< 7~pd) I v € V m} = h{V^) = 4(1^); for edges, {(I a (v), I a (u)) | 

(v, u) G E m } IAO = {(/ CT (v c ),/ ff (u c )) | (v, u) G E m } = {(I a {v),I a ( u)) | (v,u) G E’ M }\ for vertices, = 

K-pop and for elementary layers, {{L a - i (a )}Li) /<,(C) = { L J -= {^-»}a=i> because = L' a for all a. 

Now we need to show that I a (f) is an acceptable mapping for the isomorphism on the right-hand side of equation 
(2). Note that the definition of P p in equation (1) depends on the sets {L a }Qof elementary layers, and these sets are 
different in the two isomorphisms in the two sides of equation (2). We write this dependency explicitly, so that P p 
in the left isomorphism becomes P p ({L Q }g) and P p in the right isomorphism becomes P p *({L a - i( Q )}o)- With this 
notation, /<r(P p ({ia}o)) = P p ° ({L<r-i(a)}o)> so C £ P P {{La} o) => -Tt«) £ Pp- ({L*-pa)}o)- 

Now that we know that Iso p (M,M') C I a -i [Iso p ^ (A a (M), A a (M'))] for any aspect permutation a, we can 
use the aspect permutation cr _1 instead of a. Consequently, we can write I a -i [Iso p <r {A a (M), A a (M’))\ C 
I a -.\I a [lso (r) „-^A a -,(A a (M)), A a -x(A a (M')))]] = Iso p (M, M'). □ 


B. Proof of the Reduction Theorem 

We will need the following lemma for our proof of Theorem IV. 1. 

Lemma VIII.1. Suppose that f : M. —> Q c and. g maps permutations P p of M G M to permutations of G c G /( A4). 
In addition, we suppose that the following conditions hold: 

(1) f and g are injective; 

(2) /(M) 7 = f(M') =$■ 7 g g{P p ); 

(3) for all f G P p , we have f(M^) = f(M) 9 ^\ 

It then follows that Iso p (M, M') = p -1 (Iso(/(M), f{M'))). 

Proof of Lemma VIII. 1. Take any 7r G Iso p (M,M'). Because of condition (3), it then follows that f(M) 9 ^ = 
f(M n ) = f(M') and thus that g( tt) G Iso(/(M), f(M’)). This gives n G g~ 1 (Iso(f(M), f(M'))) and lso p (M,M') C 
g~ 1 (lso(f(M),f(M’))). Now let 7 G Iso(/(M), f(M')). Because of condition (2), 7 G g(P p ) and g~ l { 7) G P p . Using 
(3), we can then write that M' = / _1 (/(M')) = / _1 (/(M) 7 ) = /” 1 (/(M ) 9 ^ 9 ( 7 b) = f~ 1 {f[M 9 ( 7 ^)) = M 9 
Thus, g~ 1 (j) G Iso p (M, M'), which implies that g~ 1 (Iso(f(M), f(M'))) C Iso p (M, M'). Consequently, Iso p (M, M') = 
g- 1 (Iso(f(M), f(M'))). ' ' □ 

Proof of Theorem IV. 1. We now use Lemma VIII.1 to prove Theorem IV.1. We prove each of the three conditions for 
f p and g p that we need to apply Lemma VIII. 1. 

We begin by proving condition (1). 
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First, we show that g p is injective. Take any £, Cf £ P p such that g p {C) = g P {C)- For any a p, it follows by 
definition of P p that f a = 1 = Q' a , where 1 L a is an identity permutation over the set L a . For a £ p, the definition 

of g p guarantees that lA = v 9p A = v 9v ^ ) = tA for all v £ L a . That is, C = C,', so 9p is injective. 

We now prove that f p is injective. Take any M,M' £ M such that f p (M) = f p {M'). It follows that Vm U Vq = 
V' M U Vq, E m U Eq = E' m U E'q, and p U L Pi x • • • x L Pm = p U L' p x • ■ ■ x Lf . Because we assumed that there are no 
shared labels of vertices or elementary layers (and that tuples of elementary layers and vertices are not in the vertex 
set or in any elementary layer set), it follows that Vm = V' M and Em = E' M . Because M , M' £ M., it is also true that 
L = L' and V = V'. Thus, M = M' and f p is injective. 

We now prove condition (2). 

Consider an arbitrary 7 £ Sy fp(M) such that / p (M ) 7 = f p {M'). We want to construct £, £ P p so that g p {C) = 7 - 
For any v a £ L a , we let v% a = u 7 if a £ p and v£ a = v a if a £ p. The C defined in this way is in P p because 
permutations for a £ p are identity permutations and v a £ L a yields 7r(v a ) = a and thus vf £ L a . We now have 
by definition that v 9p ^ = v£° = vf for v a £ V 0 and vAO = for v £ Vm- If we assume that A 7 ^ v 7 for 
v £ Vm, then there exists an a such that v^ a 7 ^ (v 7 ) 0 . We know that = v a = (v 7 ) a for a £ p because of the 
coloring: (v Pl ,..., v Pm ) = 7 r(v) = 7 r( 7 _ 1 (v 7 )) = 7 r 7 (v 7 ) = 7 r'(v 7 ) = [(v 7 )p 1 ,..., (v 7 )p m ]. That is, it must be true 
that v^ a 7 ^ (v 7 ) a for a £ p. Because M' is constructed using the function f p , we know that {u, v 7 ) £ E’ 0 guarantees 
that there exists a b £ p such that u = (v 7 )f,. However, (u a ,v ) 7 = (u 7 ,v 7 ) = (iA,v 7 ) £ Eq = E’ 0 . Thus, there is a 
b £ p so that v^ a = (v 7 )^, and it thus follows that a 7 ^ b. This is a contradiction, because L a fl Lb = 0, and it thus 
must be true that v^ a = (v 7 ) a for all v £ Vm- 

We now prove condition (3). 

From a direct calculation, we verify that for all £ £ P p , we have f p (M**) = f p {M) 9p ^'. 

For vertices, we write V^ = V 9 f ^ and \J aep L^ a = \J a€p La p ^ = Vq A ^ . Combining these two equations yields 
v m C. * * * G, * * * * * AC) U V 9p{ ° = {V M U V 0 ) 9p(O = Fj p(C) . 

For edges, we write E^ = E 9 ^^ because E M c Vm x Vm, and it is also true that Eq = {(tA,A) | v £ V M ,a £ 
p} = {(u a ,v) | v £ Vm,ol £ p\ 9p ^ = Eq p(( *\ Combining the two equations yields Ej^ U E$ = E 9 ^^ U Fq p ^ ) = 
(£mU£ 0 ) 9pK) =E 9 g p{C) . 

For the color set C, the permutation Q £ P p does not change anything because it only permutes the aspects in p. 
Additionally, the permutation g p {Cf) of the vertex-colored graph does not change any vertex colors by definition. 

The color map 7 tAC)( v ) = 7 r([< 7 p (£)] - 1 (v)) = 7 r(C _ 1 (v)) = (v Pl .. -V Pm ) = 7 r(v) if v g = v £ Vm, where the third 
equality is true because C _ 1 (v) £ V M - Similarly, tt 9 ”^ (v a ) = 7r(bp(0] _ 1 (^a)) = ACH^a)) = a = ir{v a ) if v a £ L a , 
where the third equality is true because CaiVa) £ L a - □ 


C. Proof of Corollaries 

Proof of Corollary IV. 2. These results follow immediately from Theorem IV. 1. 

(1) : M = p M' <{=> Iso p (M, M') 7^ 0 ^ g~ l [lso{f p {M), f p {M'))\ ^ 0 <{=*► f p (M) S f p {M'). 

(2) : Let C be the complete invariant of = for vertex-colored graphs. That is, C{G) = C{G') G = G ', where 

G, G' £ Qc- From this invariance and (1), it follows that C{f{M )) = C{f{M')) ^==> /(M) = f(M') M =* M'. 

(3) : To obtain this result, we let M' = M in Theorem IV.1. □ 

Proof of Corollary IV. 3. The number of vertices in f p {M ) (see Definition IV. 1) is \Vm\ + \^a\, the number of 

edges is \Em\ + |Vm||p|) and the number of colors can be limited to the number of vertices. In the function f p , 

constructing each vertex, edge, or vertex color consists of copying it directly from the multilayer network or doing 

several operations of checking if an element belongs to a set that grows polynomially with the size of M. Thus, one 

can use point (1) in Corollary IV.2 to create a reduction that is polynomial in time (and linear in space) from MGI p 
to the vertex-colored graph isomorphism problem, which is known to be in GI. One can reduce in polynomial time 
any problem in GI to MGI p by mapping the two graphs to the following multilayer networks. Choose a £ p and 
use the set of vertices in the graph as a set of elementary layers in the aspect a. For the aspects b 7 ^ a, add a single 
layer 4 to the remaining elementary layer sets. For each vertex u £ V in the graph, create a single vertex-layer v 
such that v a = u and Vb = lb- (In other words, create a vertex v = (iq,..., v a ,..., Vd) = {h,. ■ ■, u, ..., Id)-) For each 
edge ( u , w) £ E in the graph, add an edge [(Zi ,..., u,..., Id), (h, ■ ■ ■, w ,..., Id)] to the multilayer network. The two 
multilayer networks are isomorphic according to exactly when the two graphs are isomorphic. □ 


